1 



O STAR X-RAY LINE PROFILES EXPLAINED BY RADIATION TRANSFER IN INHOMOGENEOUS 

STELLAR WIND 



in 
o 
o 

(N 
> 

o 



> 

0^ 



o 



Oh' 

6 



X 



L. M. Oskinova, A. Feldmeier, and W.-R. Hamann 

Universitat Potsdam, Astrophysik, 14469 Potsdam, Germany 



ABSTRACT 



It is commonly adopted that X-rays from O stars are pro- 
duced deep inside the stellar wind, and transported out- 
wards through the bulk of the expanding matter which 
attenuates the radiation and affects the shape of emission 
line profiles. The ability of Chandra and XMM-Newton 
to resolve these lines spectroscopically provided a strin- 
gent test for the theory of X-ray production. It turned out 
that none of the existing models was able to reproduce 
the observations consistently. The major caveat of these 
models was the underlying assumption of a smooth stellar 
wind. Motivated by the various observational evidence 
that the stellar winds are in fact structured, we present 
a 2-D model of a stochastic, inhomogeneous wind. The 
X-ray radiative transfer is derived for such media. It is 
shown that profiles from a clumped wind differ drasti- 
cally from those predicted by conventional homogeneous 
models. We review the up-to-date observations of X-ray 
hne profiles from stellar winds and present line fits ob- 
tained from the inhomogeneous wind model. The neces- 
sity to account for inhomogeneities in calculating the X- 
ray transport in massive star winds, including for HMXB 
is highlighted. 
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1. CLUMPED WINDS FROM O-TYPE STARS 



Young and massive O-type stars possess strong stellar 
winds. The winds are fast, with typical velocities up 
to 2500 km/s, and dense, with mass-loss rates M ^ 
lO~^M0/yr. The driving mechanism for the mass-loss 
from OB stars has been identified with radiation pres- 
sure on spectr al lines . A cor responding theory was devel- 
oped by Castor et alJ ( Il975h . After consequent improve- 
ments, this theory correctly pre dicts the observed mass - 
loss rates and wind velocities ("Pauldrach et al.', '198^. 
However, it was pointed out early ( Lucv & Solomon, 
M970), and later further investigated dOwocki & Rvbickil 
fl984) . that the stationary solution for a line-driven wind 



is unstable; small perturbation grow quickly and re- 
sult in strong shocks. The most detailed hydrodynamic 
modeling of t he line driven instability was presented by 
iFeldmeier et al.. C1997.) . These calculations show that the 
winds are strongly inhomogeneous with large density, 
velocity and temperature variations. The density inho- 
mogeneityj^s^commonly referred to as wind clumping. 
iRunacres & Owoc ki (2005) studied the 1-D evolution of 
instability-generated structures in the winds. They have 
demonstrated that the winds are essentially clumped out 
to distances of 1000 i?*. 

The theoretical predictio n of the wind clumping is con- 
firmed observationally. iBouret et alJ ll2005) conducted 
a quantitative analysis of the far-ultraviolet spectrum of 
two Galactic O stars using the last generation NLTE stel- 
lar atmosphere codes. Their study provided a strong ev- 
idence of wind clumping in O stars. It was pointed out 
that accounting for wind clumping is essential when de- 
termining the wind properties of O stars. 

One of the most significant consequences of the clump- 
ing is the reduction of empi rically derived mass-lo s s rates 
by at least a factor of 3 ( Ham ann & Koesterkei Il998t 
Bouretetal., 2005). The stellar mass-loss rate deter- 
mines the wind density, and therefore the wind opacity. 
Studying the effects of absorption in the observed X-ray 
spectra allows to probe the wind opacity. Analysis of 
the X-ray emission from individual stellar winds in the O 
stars (5 Ori (Miller et al., 2002) and C Pup (Ki-amer et al.^ 
l2003h have shown that the attenuation by the stellar wind 
is significantly smaller than expected from standard ho- 
mogeneous wind models. 

Similar conclusions are reached from the analysis of X- 
rays from colliding wind binaries (CWB). Such systems 
consist of two massive early type stars, where each com- 
ponent has a strong wind. The copious X-ray emission 
is produced in the zone where the stellar winds collide. 
At certain orbital phases the X-rays formed in the col- 
liding wind zone travel towards an observer through the 
bulk of the stellar wind of one companion. Deriving the 
absorbing column density from X-ray spectroscopy con- 
strains the wind density. Mass-loss rates can be inferred 
and compared with the models. An analysis of XMM- 
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Figure 1. Comparison between example line profiles 
emerging from a clumped and a smooth wind of the same 
mass-loss rate, velocity law ((3 = 1.5), and equal X-ray 
emissivity. Frequency is measured relative to the line cen- 
ter and in Doppler units referring to the terminal wind 
velocity Voo- The average separation between clumps h 
is indicated. 



Newton observations of t he ma ssive binary 7^ Vel was 
presented by Schild et a\\ ( l2004ft . They deduced the col- 
umn density of absorbing material from fits of the X- 
ray spectra and showed that the observed attenuation of 
X-rays is much weaker than expected from smooth stel- 
lar wind models. To reconcile theory with observations, 
ISchild et al. (2004) suggest that the ratio of clump den- 
sity to smooth wind density has to be high, and the vol- 
ume filling factor of the clumps has to be small. Signif- 
icantly, similar conclusions were r eached from Chandr a 
and RXTE observations of WR 140 (iPollock et al.Ll200l . 
Alike 7^ Vel, the column density expected from the stel- 
lar atmosphere models that account for the clumping in 
first approximation only, is a factor of four higher than the 
column density inferred from the X-ray spectrum analy- 
sis. 

Spectacular evidence of wind clumping comes from X- 
ray spectroscopy of high-mass X-ray binaries (HMXB). 
In some of these systems a neutron star (NS) is on a 
close orbit deeply inside the stellar wind of the mas- 
sive star The X-ray emission with a power law spec- 
trum is produced as a result of Bondi-Hoyle accretion of 
the stellar wind onto the NS. These X-rays photoionize 
the stellar wind. The resulting X-ray spectrum shows a 
large variety of emission features, including lines from 
H-like and He-like ion s and a number of fluorescent emis- 
sion lines. ISako et all J2003I) reviewed spectroscopic re- 
sults obtained by X-ray observatories for several wind- 
fed HMXBs. They conclude that the observed spec- 
tra can be explained only as originating in a clumped 
stellar wind where cool de nse blobs are embedded in 
rarefied photoionized gas. IVan derMeeretall (l2005h 
studied stochastic variability of the X-ray light curve in 
4U 1700-37 and its X-ray spectra. They shown that the 
feeding of the NS by a strongly clumped stellar wind is 
consistent with the observed temporal and spectral vari- 



ability. 

X-ray emission is also intrinsic for clumped stellar winds 
of single stars. Feldmeier et al. ( 1997) performed hydro- 
dynamic simulations of an O star wind. They showed that 
the de-shadowing instability leads to strong gas compres- 
sion resulting in dense cool shell fragments. The space 
between fragments is essentially void, but at the outer 
side of the dense shells exist extended gas reservoirs. 
Small gas cloudlets are ablated from these reservoirs, and 
accelerated to high speed by radiation pressure. Propa- 
gating through an almost perfect vacuum, they catch up 
with the next outer shell and ram into it. In this collision, 
the gas parcels are heated and emit thermal X-rays. The 
X-ray emission ceases when the wind reaches its termi- 
nal speed. In contrast, the cool fragments are maintained 
out to large distances, despite their expansion due to the 
internal pressure. Thus there are two disjunctive struc- 
tural wind components - hot gas parcels that emit X-rays 
and highly compressed cool fragments that attenuate this 
radiation. 



2. TRANSFER OF X-RAYS IN A CLUMPED 
STELLAR WIND 



Wind inhomogeneity alters the radiative transfer signif- 
icantly. We have studied the effects of wind fragmen- 
tation on X-ray line formation analytically in the limit 
of infinitely many stochastically distributed fragments 
(Feldmeier et al., 200 3), and numeri c ally f or more re- 
alistic stellar winds dOskinova et all 120041) . The lat- 
ter work describes the model and the code used in the 
present paper to reproduce the observed X-ray emission 
line profiles. It is assumed that X-rays are emitted by a 
large number of stochastically distributed hot gas parcels 
within the wind acceleration zone. These gas parcels 
are permeated with the numerous aligned cool shell frag- 
ments (i.e. clumps). The fragments are also distributed 
stochastically and propagate, maintaining their solid an- 
gle, outwards with v{r). They persist till large distances 
from the stellar surface. 

Because the average stellar mass-loss rate is constant, the 
mass confined within each clump is, on average, the total 
mass of the stellar wind divided by the number of shell 
fragments. The mass-loss rate M is fixed for a specific 
star from UV/optical spectral analysis. We employ the 
average number density of the fragments, n{r), as a free 
parameter. The fragments are radially compressed and 
therefore aligned. Thus the optical depth, Tsh{E, i?sh), 
of each individual shell fragment depends on its mass, 
the mass absorption coefficient at energy E, and the in- 
tersection angle between the shell and the line of sight 
(i9sh). 

The fate of each X-ray photon emitted deep inside the 
wind and propagating outwards depends on whether it es- 
capes from the wind without encountering any shell frag- 
ment, or it accidentally hits a fragment. The chance to hit 
a shell is smaller when the average separation between 
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Figure 2. De-reddened HETGS spectra of prominent O stars. All stars except Q Ori are evolved runaway stars, as reflected 
by the enhanced N/0 ratio. 



shells, h, is larger. The average separation h depends 
on the number density of the fragments and their cross- 
section: h = {n{r)a)~^. 

The probability that the photon that has hit a shell is ab- 
sorbed is 1 — exp(— Tsh). Therefore the effective opacity 
f in the stellar wind depends on the number density of 
shells n{r), the cross section of the shells cr, and the prob- 
ability of a photon to be absorbed in a shell. The optical 
depth along a photon path through the clumpy wind is 



n{r) a (1 



)dr 



(1) 



Importantly, in the limit of a very large number of frag- 
ments, Eq. Q rec overs the formu l a for t he optical depth 
in a smooth wind dOskinova et alil2004ft . 

Assume now that radiation of intensity /q is emitted 
somewhere in the wind. After passing through the wind 
this intensity will be reduced to 



/ = he- 



rn 



From Eqs. ([^E} one can see that for optically thick 
clumps (Tsh ^ 1) the effective optical depth f is inde- 
pendent of the energy. The opacity in the clumped wind 
is effectively grey. 

Our modeling shows that the wind fragmentation drasti- 
cally reduces the opacity of the wind. Therefore X-rays 



produced deep inside the wind can effectively escape, 
which would be totally absorbed in a homogeneous wind 
of the same mass-loss rate. A comparison between pre- 
dicted line profiles emerging from clumped and smooth 
winds of the same mass-loss rate is shown in FigQ 

The line width is determined by the velocity field of the 
hot plasma. It is commonly accepted to approximate the 
velocity field in the stellar wind by a "/3-law" of the form 



v(r) 



1 

1 - - 

r 



(3) 



where the distance r is normalized to the stellar ra- 
dius. The terminal velocity v^o is deduced by analyzing 
UV/optical spectra. Hence, the free parameter to describe 
the width of an X-ray line is (3. 

Intrinsic line intensity Jq can be calculated from hydro- 
dynamic wind models but it is not done yet. To compare 
the model and the observed line profiles we normalize the 
flux in the model line to the observed one. Thus only two 
free parameters are required to fit the observed line: the 
average separation between clumps h, which defines the 
attenuation; and j3, which defines the line broadening. 
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Figure 3. Thin lines: co-added meg±l lines observed in the spectrum of Q Ori. Thick lines: model line profiles. The 
rest wavelength corrected for the stellar radial velocity is indicated by Aq. Model line profiles are calculated using 
h \ R^, (3 = 2) and the appropriate for this star M and Voo- 



Table 1. Stellar wind parame ters from \ReDolust et al.\ 
lloM) and \Lamers ^tMltZOOlll 



Name 


Voc [km/s] 




CPup 


2250 


9.0 X IQ-s 


COri 


2100 


2.5 X 10'^ 


^Per 


2450 


1.8 X 10-6 


COph 


1550 


< 1.8 X IQ-'^ 



3. GRATING OBSERVATIONS OF O STARS 



We have selected Chandra HETGS observations of 
prominent O stars to fit the observed line profiles in the 
framework of a clumped wind model. The de-reddened 
spectra of selected stars are shown in Figure |2 We have 
chosen three stars that are most likely single, and there- 
fore their X-ray emission can be solely attributed to the 
stellar wind itself. The probability of a massive star being 
single is significantly higher among runaway stars and 
three stars in our sample are runaway stars that have not 
shown indications of binarity. The forth star, C, Ori, is a 
known binary. We include ( Ori in the sample because the 
initial interpretation of the observed emission lines from 
this star caused doubts in the validity of the shock model 
of X-ray production ( Waldron & Cassinelli, 2001). The 
wind parameters of the stars are listed in Table 1 . 

The observed spectra in the 4 — 24 A band were fitted 
using the xspec package. The spectrum of ( Pup has the 
highest S/N. A satisfactory spectral fit for this star can 
be achieved only by using a model that accounts for the 
continuum emission as well as for the lines. The apec 
model is capable to reproduce the X-ray emission lines 
in the Pup spectrum. The continuum can be fitted by 
thermal bremsstrahlung. Importantly, the temperatures 
obtained from the emission lines and from the continuum 
are identical, fcTx ~ 0.6 keV, and in good agreement with 



the predictions of the wind shock model. 

The plasma in hot star winds is not expected to be isother- 
mal. However, the strong shock condition implies that the 
maximum temperature of the shock material scales with 
the square of the preshock wind velocity. Constraining 
the temperature of hot gas in ( Oph provides an interest- 
ing test, because this star has the slowest stellar wind. It 
is expected that its X-ray plasma temperature should be 
a factor of two lower than in ( Pup. To estimate tem- 
peratures we used H-like to He-like ion line ratios of Mg 
and Si. Our estimates show that the temperature in C Oph 
is indeed lower than in ( Pup. A temperature less than 
0.4 keV is also consistent with the results of fitting the 
( Oph spectrum with xspec. 

The strongest emission lines observed in the spectra of 
our sample stars are shown in Figs. 3, 4, 5 and 6. Lines 
in all spectra are broad, blueshifted and symmetric. The 
Nex line is a dou blet (An,.no = 12.132, 12. 138A). 
iKramer et alJ J2003h provided a detailed analysis of emis- 
sion lines in ( Pup using the smooth wind model. They 
concluded that the lines are broadened to the terminal 
wind velocity and have similar shapes across the spec- 
trum. 



4. MODEL FITTING OF OBSERVED LINES 



A detailed modeling of the X-ray spectra and emission 
line profiles can be done only by means of full stellar at- 
mosphere models. Conventional "standard" models, e.g. 
apec, strictly speaking are not adequate for the fast mov- 
ing stellar winds, and cannot predict expected line pro- 
files. The approach of Kramer et al. (2003) was to nor- 
malize the model line flux to the observed one. The line 
profile fitting was used to constrain the wind optical depth 
at the line wavelength. This is possible because in the 
smooth wind approximation, the shape of emission line 
is sensitive to the amount of material on the line of sight. 
The more skew ed the line is, the more opaque is the wind 
lllgnacell20()l . 
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Figure 4. Observed and model NeX line in ^ Per and 
C Oph. Model line profiles are calculated using ft, « 1 
and (3 = 3 for both stars. M and Uqo are from Table 1. 



This does not hold in a clumped wind. The shape of an 
emission line is sensitive to the spatial distribution of the 
clumps in the winds. When a large number of clumps is 
effectively black, the radiation can still escape between 
them. Thus the average separation between clumps is the 
key parameter determining how much X-ray radiation es- 
capes. When clumps are compressed in radial direction, 
the line profile is distinctly symmetric although its max- 
imum is blueshifted as can be seen in Fig. 1. When the 
separation between clumps is small, the wind is nearly 
homogeneous and emergent emission lines have skewed 
shapes characteristic for a smooth wind. However, if 
the separation between clumps is large enough, the lines 
shape becomes symmetric. Therefore, the average clump 
separation can be constrained by fitting the observed line 
profiles. Shown in Fig.|3] are the lines observed in the 
spectrum of Ori. As can be seen, the observed lines are 
broad, symmetric and blueshifted. The lines are very well 
reproduced by our fragmented wind model. Line profiles 
resolved by Chandra in the spectra of ^ Per and C, Oph 
and model fits for these lines are shown in Figs. |4] and 
|5] All model fits are produced assuming that the average 
separation between the shell fragments is of the order of 
the stellar radius. 

The distance from the stellar core where X-ray emis- 
sion originates can be estimated from line ratios of He- 
like ions. The strong UV field of an O star causes ra- 
diative de-population of metastable level, weakening 
forbidden (F) line in favor of the recombination (I) line 
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Figure 5. Observed and modeled FeXNW line in ^ Per and 
Q Oph. Model parameters the same as in Fig. 4. 



dPorauet et alil200ll) . Therefore ratio F/I is a diagnostic 
of UV field. The UV field dilutes with radius, thus the ra- 
tio F/I provides information about the distance from the 
stellar photosphere. The values of F/I obtained for differ- 
ent ions are similar For instance, our analysis of Mg XI 
and Si XIII in C, Oph indicates that the plasma is located 
between 1.8 and 6 i?* from the stellar core. 

When the location of the hot plasma is constrained and 
the terminal velocity of the wind is known, parameter /3 
can be varied to adjust the model line width. As can be 
seen in Figs. 0] and |5] the width of lines differs signifi- 
cantly between stars. Lines observed in ^ Per, the star 
with the fastest wind, are much broader than the lines in 
C Oph, where the stellar wind is slowest. 

The model lines shown in Figs.|3]Sand|5lrequired rather 
large values of /3 > 2, indicating slow wind acceleration. 
This (3 is much higher than the usual /? w 0.8 inferred 
from analyses of UV/optical spectra and predicted from 
theoretical wind models (Repolust et al., 2004). The lines 
in the C, Pup spectrum, fitted with smaller values of (3 = 
1.5, are shown in Fig.|6l 



5. WIND CLUMPING IN HMXBS 



The effects of wind clumping are important for the inter- 
pretation of the X-ray spectra from the wind-fed HMXBs. 
Spectra obtained at orbital phases near X-ray eclipse are 
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Figure 6. Observed and model Ne X line in C Pup. Model 
line profile is calculated using « 0.2 _R», (3 — 1.5 and 
the appropiate for this star M and v^. 

dominated by emission reprocessed in the stellar wind. 
The X-rays from the accretion onto the compact object 
are scattered in the wind. The Compton optical depth 
scales with density, and the effects of clumping should 
be included in the modeling of the scattering in the wind. 
In addition, an emission line spectrum is emerging from 
the photoionizaed stellar wind near the compact object. 
On their way to the observer. X-rays are attenuated in the 
stellar wind. Clumping reduces the absorption and makes 
it less wavelength dependent. 

Clumping also affects the size of the photoionized region. 
Its radius is determined by the balance of ionization and 
recombination. The ionization rate scales with density, 
but the recombination rate scales with density squared. 
Therefore, if the wind consists of numerous but optically 
thin clumps, the radius of the photoinized region would 
be reduced compared to a homogeneous wind. If, how- 
ever, the wind is compressed in a relatively small number 
of optically thick clumps, then the photionized region be- 
comes more extended compared to the the smooth wind 
case (but with recombined clumps embedded). Hence, 
the radius of the photoionized region is determined by 
the clump filling factor and the density of the clumps. 



6. CONCLUSIONS 



There are firm observational evidences of stellar winds 
being strongly inhomogeneous. 

Clumping strongly reduces the opacity of the stellar wind 
for the X-ray emisison. 

Absorption in a clumped wind is effectively independent 
of wavelength. 



O star X-ray line profiles are explained by radiation trans- 
fer in an inhomogeneous stellar wind. 
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The line profiles observed in X-ray spectra of O stars are 
broad, symmetric and blueshifted. 



